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SYSTEMATIC INVESTIGATIONS OF THE INFLUENCE OF 
THE SHAPE OF THE PROFILE UPON THE POSITION 
OF THE TRANSITION POINT* 

By K. Bussmann and A. Ulrich 


The position of the beginning of transition 
laminar/ turbulent as a function of the thickness and the 
camber of the profile at various Reynolds numbers and 
lift coefficients was investigated for a series of 
Joukowsky profiles. The calculation of the boundary 
layer was carried out according to the Pohlhausen 
method which may be continued by a simplified- stability 
calculation according to H. Schlichting (4). A list 
of tables is given which permits the reading off of 
the position of the transition point on suction and 
pressure side for each Joukowsky profile. 


OUTLINE - 


I . Statement of the problem 

II. Extent of the investigation 

III. The calculation of the potential velocity and the 

practical application of the boundary layer and 
' ' stability calculations.: 

(a) Potential flow 

(b) Boundary layer and stability calculation 


*"Sy3tematische Untersuchungen uber den Einfluss 
der Profilform auf die Lage des Umschlagspunktes ." 
Zentrale fUr wissenschaf tliches Berichtswesen der 
.Luf tfahr t for s chung des Generalluf tzeugmeisters (ZVB) 
Berlin— Adler shof, Technische Berichte und Vorabdrucke 
aus Jahrbuch 19^3 der deutschen Luf tfahr tf or s chung, 

Band 10(1943), Heft 9, Sept. 15, 1943, IA 010, pp. 1-19. 
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IV. Results: .. •• • 

(a) Influence of the c a — value and of the 
Reynolds number . 

(b| Influence of the camber of the profile f/t 

(c) Influence of the thickness of the 

profile d/ t ■ 

(d) List of tables for the separation and 
instability points for all Joukowsky profiles 

(e) Mean value of the laminar-flow distance of 
suction and pressure side for all Joukowsky 
profiles 

V . Stimmary - 

VI. References 


SYMBOLS 


x,y rectangular coordinates in the plane 

s profile contour length starting from 

the nose of the profile 


t 


wing chord 


^t» 


U m (s) 

6 P 4 


length of the profile contour from 

nose to trailing edge (different for 
pressure and suction side) ' 

velocity of. incoming flow 

potential velocity at the profile 

boundary layer thickness according to 
Pohlhausen P4 


6 * 


z 4 

*P4 



— * ' Ji- 
" z 4 U 0 W 


displacement thickness of the boundary 
layer 

nondimens.ional boundary layer thickness 

form parameter of the boundary-layer 
profiles according to Pohlhausen P4 
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s crit 
S AP 6 


: .form parameter according to. 

;• -Pohlhausen P 6 

universal' functions of the boundary- 
layer calculation 

position of the instability point, 
•measured along the contour of the 
: nose of the profile 

position of the separation point 
. according to P 6 method . • 


I. STATEMENT OF THE PROBLEM 1 


The position of the transition point laminar/ turbulent 
in the frictional boundary layer is of decisive importance' 
for the problem of the theoretical calculation of the 
profile drag of an airfoil since the friction drag 
depends on it to a high degree. The position of the 
transition point on the airfoil is largely dependent 
on the pressure distribution along the contour of- the 
profile and, therefore, on the shape of the airfoil 
section and on the lift coefficient. A way of theoretical 
calculation of the start of transition (instability 
point) ^ that is, the point downstream from which the 
boundary layer is unstable, was recently indicated 
by H. Schlichting (1,3,4) and J. Pretsch ( 2 ). 

According to present conceptions the turbulence 
observed. in tests develops from an unstable condition 
by a mechanism of excitation as yet little known; ' ■ 

therefore, the experimental transition point is always 
to be expected a little further back than the theoretical 
instability point. 

Knowledge of the. theoretical instability point is, 
nevertheless, important for the research, on profiles, 
in particular for the drag problem. Recently a report 


An extract of this report was given in a lecture 
of the first— named author at the Lilienthal meeting 
for the discussion of boundary— layer problems in 
GOttingen on October 28 and 29 , 1941, 
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was made. about airfoil sections whi-ch, due to a position 
very far back of the instability and .transition point, 
have surprisingly small drag coefficients (laminar 
profiles). Thus . far no systematical investigations of. 
the influence of the shape of the profile upon the 
position of the transition point have been made either 
experimentally or theoretically. The following calcu- 
lation of the., theoretical instability point is, 
therefore, given for the first time in a sufficiently 
large range of c a — values and Reynolds numbers to 
achieve a greater systematization of airfoil sections. 

In order to keep the extent of calculations within 
tolerable limits only the two most important profile 
parameters, thickness and camber were varied. A. rather 
convenient and accurate mode of calculation of the 
potential flow for the profiles is important for> these 
investigations and the selection of a series of 
Joukowsky profiles was, therefore, • natural . It was 
not advisable to take for instance the NACA series as 
a, basjs ; the calculation of the potential' flow for 
such profiles according to the methods at present 
available does not achieve the accuracy which is 
required here . 


II. EXTENT OF THE INVESTIGATION 


A series of ordinary Joukowsky profiles of the 
relative thicknesses d/t = 0 , 0 . 05 , 0 .10, 0 .15, 0 .00, 
0.25 and the relative cambers f/t = 0 , 0 . 02 , 6.04, 

0.08 were taken as a basis, (See fig. l.j For instance, 
J 415 stands for the Joukowsky profile of camber 
f/t = 0.04 and the thickness d/t - 0.15. The c -region 
which was examined is c„ = 0 to 1 and. the Re —number 

range Re - -^ = 10 + to 10 . The complete calculations 

were carried out only for the following profiles: 000 , 

005, 015, 025, 215, 400, 415, 425, 300, 315, and 825 . 

The results for the remaining. profiles could be obtained 
by interpolation. Thus it was possible to obtain a 
result with tolerable loss of time in spite of the very 
extensive program (four parameters); a certain amount 
of accuracy had to he neglected since the interpolation 
sometimes 'was carried out over three points . 
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III. THE CALCULATION OP THE POTENTIAL VELOCITY 
AND THE PRACTICAL APPLICATION :OF THE' BOUNDARY 
I AYER- AND STABILITY CALCULATIONS' ' 

, t 

(a) The Potential- Flow 


The calculation of the potential velocity with its 
first and second- derivatives along the profile contour 
forms the basis for a boundary layer and stability 
calculation. The potential flow about a Joukowsky 
profile is obtained by conformal mapping of the flow 
about a circular cylinder. (See fig, 2 .) ' 


A short list of the most important symbols and..' 
formulas for the profile contour and for the velocity 
distribution follows: 


z = x + iy 
S » 6 +. irj 
mapping function: 




*>= Coordinates in the complex plane 


I = 


z + 


circle > mean camber line of the prof ile 
circle K — J cambered profile 

a . radius of the unit circle in. the z— plane . - 
R radius of the qircl.e to be mapped in the z-plane 
t wing chord 


-• length of the profile contour from nose to 

trailing edge (different for suction and 
pressure, side) 

x c>f 7 0 - coordinates of the center of the circle to 
be mapped in the z-plane (circle K) 
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o, n 

<P 

-P 

a 


a 


g 


center coordinates of the:. mapping circle of the 
mean camber line of the' profile (circle K4) 

varying angular' coordinate of the conformal 
transformation' 

zero lift direction (See fig, 2.) 

angle of attack of the airfoil referred to 
the theoretical ■ chord 

geometrical angle of attack referred, to the 
bi tangent (See fig, 2.) 


Profile nose: cp = it + p 

Trailing edge: <p = .-3 


x 


o 


ss = thickness parameter? k = 1 + 


? See table I . 


71 


~ = = camber parameter 

1 


P = arc cos 




1 + e 


2 


= arc sin 


J 

e 4 


4 .. 


\ 


■ 2 
1 + € 4 


The profile parameters and P can be found 

in table 1 . 


Profile contour: 

r 

k 


x _ a 
t “ t 


1 + 6 ), cos cp - e 


:i 


i +■ |\ + 


-2\ 

2 


k + «. 


Z - a 
t " t 


k f e 4 + \j 1 + s 4 2 sln 9 ) f 1 " I 


>(l) 


J 
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Nose radius 9/ t: 

For symmetrical profiles the equation' 


P 

t = 


2e 


1 


i d 1 + 4® # 


( 2 ) 


is valid exactly. This formula may with a good approxi- 
mation also he applied to cambered profiles. The 
numerical values in table 1 show that the nose radius 
of the Joukowsky profiles is' only little larger than 
for the NACA profile family according to NACA report 460_ 
for which P/t = l.l(d/,t) . 


a) + sin (a + p)J P 1 ( cp ) ( 3 ) 


Velocity distribution: 


U. 

U 


m 


= 2 


o 


jja in ( cp 


Stagnation points: Back q> = — 0 

Front © = 7 r + p + 2a 

' - N 




( 4 ) ' 


\}(W - l ) 2 + 4k 2 (e ^ + ^ 


1 + 6 ^2 s ^ n m ^2 


<pY 
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Velocity at the trailing edge: 


lira / cos (a + p) 

cp 


o) 


“f 1 + \ S 


( 5 ) 


Arc length: 


ds 

dcp 


\ 2 


\(§) 



( 6 ) 


s as a function of cp is. to he ascertained from (6) by 
graphical integration or can be seen directly in an 
enlarged presentation of the profile contour" (t = lm) . . 

Velocity gradient: 


1 dU : 


m 


u 


0 dcp 


2 


N 


N, 


V? 


Cos (cp — a) 


1 


+ — §/2’{^ in (<P ~ «) + sin (a - p)j 


N 


A = : 2k 


+ i 


1 + sin 


cp) 


N x = (N - l) 2 + A 2 


B = 


N x N* - N 


0 - 1) 


N» + 2k 




2 ~] 
1+6. A cos o 

4 . u 




N r = 2k \J1 + «/ A sin cp + € k cos 

4 V 1 4 


’) 


( 7 ) 
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The first derivative of the potential velocity with 

d.U m 

respect to the arc length "g~ s ~ was calculated 
numerically from equations (6) and {7), and from that 
graphically the second derivative 




ds 2 * * * 


Relation between c_ and a: 

cl 


R 

c a = 8 tt -.j--sin (a. + p) 


R 

t = 


kyl + «4‘ 


3 + 2 6 -j + y 


1 


j compare tables 1 and 2. 


k + € 


(b) Boundary— Layer and Stability Calculation 


After calculation of ; the potential velocity with 
its first and second derivatives along the profile 
contour there is a boundary— layer and stability calcu- 
lation to be made for each profile. The boundary— layer 
calculation according to- P’ohlhausen (5) was based upon 
the differential equation for the boundary layer thick- 
ness in the shape indicated by Howarth (6)*. • ' , 




z 4 2 g (*)U» 



2 In the meantime a simpler form of the Pohlhausen 

equation was indicated by H. Holstein and T, Bohlen (10) 
where the momentum thickness appears as independent 

variable. For this method the second derivative U" 

is unnecessary; the integration procedure is thus 

simplified considerably. 
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Hie following - s^iibo is "stand fbr ' : 


__ u m 

U ~ tj ” > U* 

o 


.117, 


j.2 d^I.L 


TT „- dm... 

If ds • » ^ Z U ~T> 
0 o a s 


>pj, - boundary layer thickness according to Poillhaus'en 


si. 


°B4 2 Tj o 


\ 


p|,- - 2^ T J» = form parameter according to Pohlhausen 


■(^pij) and g(\ pJ ) = universal functions 




. lib. / 2 1 V 2 2 

2 ~ 51^' loirs + 120r' • + Wr 


x5 


■ XL 

b 6 0 




■A 2 


fc>yo bpu icp.44 


1 


g(\) = 


94 S 


+ 2 X . 


.2222. 


12. . 




/ 70 / -V ,-x 

O^o fW! 


J2L 


630 iop.44. 


Initial conditions s ' 

. At the stagnation point 


\ = 7.052. 


that is.. 


z, 


40 


_ 4o = 7 .052 


U» 


U» 


~ (9) 
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11 ' 


Besides , 



U" 

- 5.391 ~T 2 


( 10 ) 


The isocline method was selected for the solution of 
the differential equation. The particular advantage of 
this method is that not. only the initial value z^q is 

known, but that the initial inclination at the stag- 
nation point also can be determined. The latter 

value is obtained by exact performance of the limiting 

process lim |f in (8) (I-Iowarth(6) ) . With z^' known 
U-i 0 , ' . 

the integral curve passing through the initial value z^q 
is easily found which otherwise is not immediately 
possible" because of the singularity of the Pohlhausen 
equation at the stagnation point. 

For the profiles of the thickness d/ t = 0, that 'is, 
for the flat plate and the circular-arc profiles, the case 
where the flow does not enter abruptly (a = 0) is 
exceptional since there exists no true stagnation 
point: the velocity at the leading edge has a finite 

value different from 0. The initial value of the 
thickness of the boundary layer is here zero, that is, 
at the leading edge there is: 

' v . i - 0 ' . (n > 


Profile 

c not abrupt 
a flow entrance 

0 

0 

200 ; 

.25 


} 

400 

\ .5 

800 

1 
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The velocity near the leading edge of circular- arc 
profiles takes the same course as \f s: 


^m- ^00 + C '/s + . . . 

/ 

that is, U* becomes with -s-^ 0 infinite like 1/ \flT- 
the velocity has a perpendicular tangent' which always 
occurs when the contour- of ..the profile shows a sudden 
change in curvature as it does, here (v. Koppenfels (8).). 
zij. near the leading edge for a circular-arc profile 
behaves like z^ for the flat plate, that is, goes 

5.11 a linear relation to .. s. toward 0. 3 Taking these 
facts into consideration there results at the leading 
edge: ■ 


1 

o 



34.05 

U oo 


( 12 ) 


It has proved advantageous to calculate the line 
elements zs/ directly from the equation (8) by means 

of a plotting of the curves 


and g(* pJf ). (See 


fig. 3.) This method is superior to the calculation of 
the line elements by means of the often used nomograms 
of Mangier ( 7 ) with respect to accuracy and its equal with 
respect to loss of time. Generally it will be sufficient 
to determine the line elements for each value of the 
abscissa s/t at two ordinate values only. 


The boundary-layer calculation yields for each 
profile for a given c 0 value the nondlmens ional 

boundary layer thickness z^ and the form parameter 

as a function of the length of the arc s along the 
contour.. The distribution of velocity u(y) in the 
laminar boundary layer is then! obtained from: 


For the flat plate, z, = - 34.03 s/t 
Pohlhausen (5)). ^ 


(according to 
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The results of the boundary-layer calculation for 
the profiles J 800 and «T 025 have been plotted as 
examples in figures k _and 5 : the form parameter X ^ 

and the nondime nsional displacament thick- 


ness 



with 6-"* standing for the displace- 


ment thickness. 


The following relation exists between the displace- 
ment thickness and the boundary layer thickness according 
to Pohlhausen: • . ' ' 



The, displacement thickness of the flat plate in longi- 
tudinal’ flow ^pji = 0^ is represented graphically 

in figures k ^nd 5 for comparison* The following equations 
are valid: 



6* 5 P1l JM - 
6 Pk \ V ” 



(16) 



r 
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• The profile J 800 (fig. 4)- ..shows clearly that 
the displacement thickness for-! accelerated, flow (suction 
side) is smaller than the displacement thickness of the 
flat plate whereas it is. larger for retarded flow 
(pressure side). (Compare also fig. l6.) 


Prom the boundary- layer calculation there result 
also the laminar separation points,.. According to 
the four -term method of Ponlhausen separation occurs 
at Ap[|_ = -12, according to the slut- term method (see , 

below) at »x p £ = -10 corresponding to 


• . ..( 

Plow photographs have been taken in a Lippi sch 
smoke tunnel for a 'part of the. calculated profiles of 

tJ o t 

models of 50-centiraeter w ing chord and at Re -numbers v 

of about 2 x 10 5 . The points of separation have been 
ascertained from the flow graphs (figs. 6 to 11, appendix). 
Figure 12 shows the experimental and theoretical separation 
points for various profiles for comparison. Compare 
also table 3, The agreement is rather good. 


After Xtd| ( has been ascertained as a function of 

the length of the arc s , there results the instability 

point (s/t) from a stability calculation 

cr it# ^ 

(H. Schilchting (k) ) based, on- the six-term method of 

Pohlhausen. The P6 -method is based on a one-parameter 

group (parameter X.p6 ) ot boundary- layer profiles 

which can be represented by polynomials of the sixth degree. 
An investigation of stability was carried out for.a * 
number of these boundary '-layer profiles in (i|)j first, 

the critical Re -number of the boundary layer “TJ I 

• ■ . crit. 


as a function of Apg was obtained. The. critical 


/j/' 


Re-number of the laminar layer 1 as a function 

\ . 'crit. 


of Xp| + ’ (fig. 13) is then immediately known .also - 
because of a universal relation between Xp6 an ^ ^ pip 
indicated in (ij.). 
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Once \pj^(s/t) has been ascertained from the 
boundary-layer calculation according to Pohlhausen ' s 


method a critical Re-number : "\ 


IV m 6 


i 


may • be 


\ v /trit • 

coordinated to each point of the profile by means of 
figure 15 r> Moreover the • Re -number of the boundary 


layer 


• T r k# 

u m^ 

u 


can be calculated for each point of the' 


profile at a certain 


V 


u 6 

m 




TL 


6*’ JV 


ill J.H V . > 

— = n b- 


(17) 


The- location- of the instability point is then- given - by - 


U 5 
m 


* 


/ 


v 


=l 


'U 6 f 
m 


V v / crit* 


(18) 


‘ IV. RESULTS 

(a) influence of the c a - Value and the R 9 -number 


The -results of the stability calculation, that is, 
position of the theoretical instability point /-f^ 

' '%rlt . 

for the sample profiles J 800 and J 025 are plotted 
in figures lip and 15 against c & with the Re -number 

U G t 

as parameter and furthermore against -d— with the 

■ v 

c. g .- value ; as parameter. The' characteristic course of the 
curves, is the same for all profiles j the following 
statements are valid! the. instability point travels, 
with increasing* c a at a constant Re-number, forward 


the 
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on the- suction side, backward oh the ‘pressure side; 
the instability point- travels forward on both suction, 
and pressure sid'e with increasing Re -number ' at a'' 
fixed c a -value. y This behavior is demonstrated very 
clearly in figures l6 and 17 which represent the 
velocity distributions for the two profiles J- 800 and . 

J 025 for the various c a -values' with instability and 
separation points. One .can see in particular, that the . 

» instability* do ints of the suction side for Re-numbers 
U t “ •_ 

from ' ■ ~ 10 ^ to 10' lie near the velocity maximum; 

v , 

mostly the position of the instability point for Re = 10 
agrees well with -the location of the velocity maximum. 

The pressure side of J. 800 .in the case where the flow . 
does* not enter abruptly (c a = 1) is an exception 
among the above mentioned examples, since the flow from 
the leading edge to the center of the profile is con- 
siderably increased so that no relative velocity 
maximum exists. Measurements concerning the dependency 
of the transition point on the c & - value were taken 
by A. Silverstein and J. V. Becker (9)» These tests 
showed (as a result) the same dependency of the 
transition point upon the lift coefficient as the 
present theoretical investigations. 


(b) Influence of the Camber of the Profile 

’ The influence of the camber upon the position of 
the instability point can be described as follows § 
the instability point travels with increasing camber, 
at. constant thickness, for all c a -values and Re-numbers 
backward on the .suction : side, forward on the pressure 
side. This influence of .the camber can be understood 
from the fact that the stagnation point and therefore 
the region of the accelerated stabilizing flow travels, 
with increasing camber, backward on the suction side 
whereas because of the flow around the no3.e of the 
profile a region of considerably retarded destabilizing 
flow originates immediately behind the nose on the 
pressure side. Figure l8 represents as an example the • 
results for profiles of the thickness ■ d/'t =/0*15, 
with variable camber f/t for c = 0.25 v and again 
the Re-number as parameter. The curves for all thick- 
nesses and all c Q -values have the -same characteristics# 

cl 
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• ••■ • . (c) Influence of. the. Profile Thickness; : 

The dependency of the instability point on the 
thickness cannot be described in such general terras 
as the influence of the camber since this influence 
depends in the following way on the c a value : A 

-certain .. 11 c a no +. a b ru pfc flow entrance ” » that is, the - 
c Q value that corresponds to the not abrupt entering 

a. .. -7 * . - ' 

of the flow ( a - 0) for the circular arc profile ; ; 

with the given camber , i3 coordinated to each value of 
the cambeh f/t. The curves. ( s /t) cr it, versus d/t , 

at 'a constant f/t show on principle two -different 
types (fig, 19 ) : - ; : ' 

I. With increasing thickness, the curves (s/t) Qr ^r 

versus d/t start from a finite value and have a 
flat minimum: < 

, On the suction side for c a - c a for not abrupt 
flow .changes. 

On the pressure side for c a - c a for not abrupt 

flow changes, 

\ 

' II. The curves :(s/t) cr;Jt versus d/t rise starting 

from 0 with increasing thickness; hence, the transition 
point moves backward as follows: 

On' the suction side for c a > c a for not abrupt 
- flow changes, 

. ; On the pressure side for c a < c a for not. abrupt 
flow.;, change s . " J . ' 

1 ihh. results for the symmetrical profiles at c a - 0,25 
are represented as an example in figure 20. For the 
symmetrical profiles' c a t or no t abrupt flow changes = 
.that i3, the dependency of the instability point on the 
, thickness d/t for all c a > O' is of type II on the 
suction side, of type I. on the pressure 'side , . 

The flat minimum in curves of type I does, in some ~ 
cases, not exist at high. Re -numbers (Re =10? to 10°), 
and (s/^Jerit versus . -d/t rises' from- the finite 

value d/t = 0 « v " • • • : • • " *•%/* . 
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(d) List of Tables for tlie’ Separation and Instability 
Point in all Joukowsky Profiles 

The total result of the boundary-layer and ' 
stability calculations is represented by a graph of the. 
curves (s/t) Ap £ = const., and (s/t) crii _ . = const., respec- 
tively, in a system of axes thickness d/t - camtier f/t. : 

(See figs, 21 to 30, ) A, profile corresponds to each point 
of the plane. In particular, the : 'symmetric al profiles 
are' coordinated to the Points of the d/t -axis, the 
circular arc prof iles to the. points, of the f/t-axis, 
and the flat plate corresponds to the zero point. . : 

Lift coefficient and Re -number are considered as 

parameters.. One has therewith a catalogue of Joukowsky 

profiles that make it possible read off ,. for every 

profile in the region 0 = d/t - 0,25;. 0=f/t - 0/08, 

the position of the separation points for 0 '= c - .1 
(figs, 21 and 22) and the position of the ins tab ill ty.. r -.. 

point for.’ 0 = ca = 1 and 1C)5 = Re - 10°. Figures 2$ to 30 

represent the curves . (s/t . •*. = const, for the ' 

CPXut. 

U n t c, p. 

Reynolds numbers from Re - — ~ - - 10-' to 10° at 

the . c a -values c a - 0, 0 . 25 # 0 . 5 * and 1 for suction 

and pressure side. For instance the values indicated 

in the following table for profiles of the camber f/t = 0,02, 

and the thickness d/t = 0,10 to O.15 at Re = 10^ and 10^ 
are taken from these representations. (See page- 19 ,)' 

The most remarkable matter in this graphical repre- 
sentation is the location of the curve (s/t) A p £ = 0, and " 

(s;/t) cr ^^ > =.0, respectively, at the various. c a - values-, . 

The position of this zero curve in the catalogue., for the 
instability points will be discussed; the same is valid 
for the separation points* ^A^crit. ~ ® can ° n l7 , 

appear for the flat plate and the circular-arc profiles 
on the suction side for c a > c a for not abrupt flow 

ohanges, on the pressure side for c a < o a for not 
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. The curve (s/t) =0 always 

abrupt flow changes . crit. 

coincides with the f/t-axisj it forms a part of 

the f/t-axis which- is determined by the actual c a -value. 

Therefore no point (s/t) ' * =0 exists on the suction 

CI*1 1 e 

side xor c a ~ 0 oince f Qr no ,_ al3rU pt fi ow changes ^ 

for all circular-arc profiles# For the pressure side, 
on the other hand, (s/t) : =0 bn the whole f/.t-axis# 

crit. 

There follows in the same vj ay for c a = 0.2J 

that (s/t) .. =0 for 0 = f/t < 0.02 on the suction 

7 crit. 

side and for f/t > 0.02 on the pressure side. Pressure 
and suction side, therefore, always complement each 
other. The point which corresponds to the circular-arc 

profile with c- a - ° a f or no t abrupt flow changes 

(for instance J it00 at c & = 0.5, compare figs. 2J to JO), 

that is, the end point of the distance (s/t) =0 

is a singular point in. the following sensei The point 
itself assumes a certain value (s/t) (different 

crit. 

for pressure and suction side), but an infinite number 
of curves (s/t) =- const, which are crowding 

crit. " . 

together asymptotically toward (.s/t) , . <= 0 run 

into it. It is true, these relations for the very thin 
profiles give only qualitative results from the present • 
investigations. An additional series of thin profiles 
would have to be investigated in order to make more 
accurate statements possible. -However, only profiles 
with thicknesses d/t > O.OJ which can -be analyzed 
quantitatively, are of practical interest. 

For c =0, (s/t)' . is the same on suction and 

a * crit. 

pressure side for the symmetrical profiles. Therefore 
the curves (s/t) = const, for suction and pressure 

' crit., 

side would adjoin at c a = O in a joint representation 
of the suction and pressure side where for the pressure 
side the measure of the camber is directed downward. 

For values c Q / 0 also the curves (s/t) = const. 

1:1 , crit. 
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have continuations which correspond to the 
curves (s/t) = const* for the pressure and 

suction side, respectively, at the appertaining 
with inverted sign. ", ^ 


,- value 


(e) Mean Value of the Laminar- Plow Distance on 

f ... 

Suction and Pressure Side for all Jcukowsky Profiles 

In view of the development of laminar profiles 
tne mean value of the laminar- flow di stance on suction 
and pressure side is interesting, Figures 31 and 32 " 
show the curves mean value (s/t) ”= const, in 

C ^ 1 t 4 

the d/t-, f/t -plane for various lift coefficients and 

the He -numbers — ~ =10° and 10? jV '=0.3 
> ’ u ■ _L crit . y 

l^crit . suction + S crit. pressure )| . L * ene ?ally the 

following conclusions are valid: The profiles with the 

smallest mean value (s/t) ; for a certain c -value 

r 1 u • 8. 

lie near the circular- arc profile to which this value 
is coordinated as c 

a for not abrupt flow changes* 

This profile will be for c & = 0 the flat plate, 

for c a = 0.23 the profile J 200 , for c a = O.p the 

profile J t .00 and finally for c a = 1 the profile J 800 


There seems to be an exceptional case. at he 


Vt 

v 


= 10 ( 


and c a - 0,5 (fig. 31) which can be explained as 
follows: The circular-arc profile for which at the 

considered c a - value the flow enters "not abruptly" 

(for instance J 4OO at c & = 0.5) is a singular point 
in the f/t-, d/t-diagram 7 Approaching this profile 
on the f/t- axis from two different sides. one obtains 
two dinerent limit values (s/t)^^^^ , since once only 

the suction side and once only the pressure side 
contributes to the mean value. Only for the singular 
point itself suction and pressure side. contribute so 
that ohis profile has a higher (s/t) than the 

crit. 
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profiles on the f/t-axis_ : hear it. If one now considers 
the curves mean value (s/t)^.^ - const, for values 

higher than the two. limit values (range I). These 
curves enclose the singular point and end at two points 
on the f/t-axis. The remaining smaller mean 

values (s7t) generally cover only a small region 

cr.it . 

near the singular point (range ll) where, with the 
present investigations as a basis, more accurate 
statements are not possible. Only for the, 

U n t f 

esse — — - = 10° and c = .0,p the range II comprises 
v a 


all profiles of the series considered here since on the 
pressure side the profile J kOO at c = Q.p and 


He = 10 b has no' transition point 


<s/t) crit. = 1 


and 


point f/t = 0,0k obtains a high mean 

> 0.5. For this' case there are closed 
crit. . 

curves (s/t) = const, and there exists a profile (J 115) 


therefore. th 
value (s/t) 


with the smallest mean value 
Re = 10 6 , 


Went. = °’ 155 at 


Moreover, the following results are obtained from 

figures pi and 7 /l\ All Joukowsky profiles have small 

mean values (s/t) ; for instance, the mean values 

' crit. 

for practically Important profiles with the camber 
f/t = 0.02 and the thickness d/t = 0.10 to 0.20 at 

Re-numbers of 10° to 10^ are between 0.08 and 0.2. 

These mean values are only to a small degree dependent 
on the lift coefficient; for Instance, the me,an values 
for the profile J 215' at Re = 10° and at lift coef- 
ficients c = 0 to 1 are between 0.l6 and 0.175* 
a 


V. SUMMARY 


A series of Joukowsky profiles with thick- 
nesses d/t = 0 . to 0;. 25 and cambers f/t =. 0 to 0.08 
was investigated with respec.t to the position of the 
instability point for various lift coefficients and 
lie-numbers. The following result was obtained j v/ith. 
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increasing Re-number, the instability point moves 
forward on suction and pressure side; with increasing 
c a - value- it moves forward on the suction side, backward 
on the pressure side, . The position of the' instability 
point as a function o.f thickness and' camber of the 
profile is represented in .the shape of a graphical 
list' of tables which permits the- reading off of the 
position of the instability, point on suction and pressure 
side as well as of the mean value of' .the laminar-flow 
distance on suction and pressure side for- each profile 
of the series , 

Translated by Mary L, Mahler * ■ • . 

National Advisory - s 

Committee for Aeronautics 
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TABLE 2 

THEORETICAL ANGLE OP ATTACK a (DEGREE) 

c a 

Profile 
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000 200 ICO 000 



Figure 1.- Joukowsky profiles: thickness d/t = 0 to 0.25; camber 
f/t = 0 to 0.08. Profile number: for instance, J 415 stands for 
the Joukowsky profile with f/t = 0.04 and d/t = 0.15. 



Figure 2.- Explanatory sketch to the Joukowsky transformation (schematic) 
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Figure 3. - Auxiliary function f ( *■ p4 ) . and g(^ p4 ) for the boundary-layer 

calculation. 



Figure 4 - Profile J 800, boundary -layer calculation: form parameter \ p4> 


and displacement thickness 




t 
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Fierure 5.- Profile J 025, boundary -layer calcualtion: form parameter \ 

5 li 


and displacement thickness 
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a* 

C a tlieor. 

—4,6 

0 

—2,3 

0,25 

0 

0,50*) 

2,3 

0,75 

4,6 

1,00 


a* 

Ca theor. 

—9,2 

0 

—6,9 

0,25 

—4,6 

0,50 

— 2,3 

0,75 

0 

1,00*) 


Figure 6. Profile J 400, smoke tunnel 
photographs (Re = 2. 10*). 

*) not abrupt entering. 


Figure 7. Profile J 800, smoke tunnel 

photographs (Re = 2. 10^). Laminar 
separation points see Fig. 12. 

*) not abrupt entering. 
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a * 

Ca tlwor. 

0 

0 

2,2 

0,25 

4,4 

0,50 

6,6 

0,75 

8,8 

1,00 


Figure 8. Profile J 005, smoke tunnel 
photographs (Re = 2.10^). 
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a # 

Ca tbeor 

0 

0 

1,9 

0,25 

3,8 

0,50 

5,7 

0,75 

7,65 

1,00 


Figure 9. Profile J 025, smoke tunnel 

photographs (Re = 2.10 5 ). Laminar 
separation points see Fig. 12. 
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a* 

V 

Ca tbeor. 

— 4,6 

—8,1 

0 

— 2,5 

— 6,3 

0,25 

— 0‘,5 

— 4,2 

0,50 

1,45 

— 2,2 

0,75 

3,5 

— 0,2 

1,00 


Figure 10. Profile J 415, smoke tunnel 

photographs (Re =-2.10^). Definitions 
of oc and oc g see Fig. 2. 


a * 

* g ° 

Ca tbeor. 

— 9,2 

— 11,3 

0 

— 7,2 

— 9,3 

0,25 

— 5,7 

— 7,25 

0,50 

— 3,2 

— 5,25 

0,75 

— 1,1 

— 3,25 

1,00 


Figure 11. Profile J 815, smoke tunnel 

photographs (Re = 2.10^). Laminar 
separation points see Fig. ]2. 



NACA TM No. 1185 


39 


tei _ . . ., Theoretical 

(J/ — Suction side 

(D) — pressure side — Q— Experimental 



Figure 12.- Laminar separation points A P g versus c & , comparison of test 
and calculation for the profiles J 800, J 025, and J 815. 
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and the form parameter X 


P4 * 
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i’ 




Figure 15.- Profile J 025: 


Result of the stability calculation, ^ 


Uot 

u 


and c 


a* 


versus 
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Figure 16. 


\ 



- Profile J 800: Velocity distribution with instability and separation 


points at various 



and c a - values. 



BV K"* 

is tability points 



Suction side 
pressure side 


1 



co-0 

—^^1025 
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v« 


to 
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Figure 21.- Position of the laminar separation point (s/t) . _ _ as a function 

AP 6 

of the thickness of the profile d/t and the camber of the profile f/t; suction 
side. 




Figure 22.- Position of the laminar separation point (s/t) Ap g as a function 

of the thickness of the profile d/t and the camber of the profile f/t; pressure 
side. 
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Figure 23.- Position of the instability point ( s /t) cr p- as a function of the 
thickness of the profile d/t and the camber of the profile f/t; suction 
side; Re = 10 5 . 



Figure 24.- Position of the instability point (s/t) cr j t as a function of the 
thickness- of the profile d/t and the camber of the profile f/t; pressure 


s: 
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Figure 25.- Position of the instability point (s/t) crit as a function of the 
thickness of the profile d/t and the camber of the profile f/t; suction 
side; Re = 10 ^ 



Figure 26.- Position of the instability point (s/t) crit as a function of the 
thickness of the profile d/t and the camber of the profile f/t; pressure 
10 6 . 
















